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We investigate a counterintuitive geometric interaction between defects and curvature in thin 
layers of superfluids, superconductors and liquid crystals deposited on curved surfaces. Each defect 
feels a geometric potential whose functional form is determined only by the shape of the surface, 
but whose sign and strength depend on the transformation properties of the order parameter. For 
superfluids and superconductors, the strength of this interaction is proportional to the square of the 
charge and causes all defects to be repelled (attracted) by regions of positive (negative) Gaussian 
curvature. For liquid crystals in the one elastic constant approximation, charges between and 47r 
are attracted by regions of positive curvature while all other charges are repelled. 

PACS numbers: Valid PACS appear here 



The physics of topological defects on curved surfaces 
plays an increasingly significant role in the engineering 
of devices based on coated interfaces. 0,l3il3- Defects 
also affect the mechanical properties of some biological 
systems, such as spherical viruses, whose shape is de- 
pendent on the presence of disclinations in their protein 
shell 3. Furthermore, the effects induced by a curved 
substrate on the distribution of defects are not fully un- 
derstood even in well studied systems such as thin super- 
fluid or superconducting films. In this paper, we study 
simple continuum generalizations of the plane XY model 
to frozen surfaces of varying curvature to gain a broad 
understanding of the interaction between topological de- 
fects and curvature. 

The XY model is a simple setting in which particle-like 
objects emerge from a more fundamental theory. The ba- 
sic degree of freedom is an angle- valued function on the 
plane whose values vary from to 27r. These angles could 
represent the orientations of interacting arrows. The in- 
teraction, which tends to align neighboring arrows, re- 
sults from the continuum free energy T given by 



T^^ J£n{V9{n)f , 



(1) 



where the set of coordinates u = {x,y) label points on 
the plane. Despite its simplicity, this model captures the 
main properties of vortices in layers of superfluid *He or 
thin superconducting films when the field (u ) is iden- 
tified with the phase of the collective wave function. In 
addition, the elastic energy of Eq.lQJ correctly describes 
liquid crystalline phases for which the bond angle, ^ (u ), 
has periodicity — with p > 3. For a solution of nemati- 
gens {p = 2) and tilted molecules in a a Langmuir film 
{p = \) two different elastic constants are necessary to 
account for bend and splay deformations 0, but these 
are renormalized to the same value at finite tempera- 
tures Q. Besides its experimental significance, the XY 
model is the cornerstone of our conceptual understand- 
ing of topological defects, singular configurations of the 
field 6* (u). 

Like particles, defects have charges and a characteristic 
Coulomb-like interaction. The charge q, a multiple of 



— , can be defined by the amount increases along a 
path enclosing the defect's core. The force between two 
defects located at positions u^ and u^ is determined by 
the energy Eint stored in the 9 field: 



Eint = K Qi QjU {u.i,Uj) 



(2) 



where the inter-particle potential U{ui,Uj) is propor- 
tional to the logarithm of the distance in the plane. 

On a fiat surface, thin layers of superfluids, supercon- 
ductors and liquid crystals can all be analyzed within the 
framework of Eq.JQ) 7]. However, there is a crucial dif- 
ference between, say, the phase of the superfluid order 
parameter and the angle that describes the local orienta- 
tion of liquid crystal molecules. The former transforms 
like a scalar since the quantum mechanical phase does 
not change when the system is rotated, while the lat- 
ter represents a vector aligned to the local direction of 
the molecules. Thus, a common boundary condition for 
a liquid crystal is for the director to be tangent to the 
boundary of the substrate. By contrast, no such con- 
straint exists for a *He film because its wave function is 
defined in a different space from the one in which the 
superfluid is conflned. This distinction is crucial on a 
curved surface. In the ground state of a ^He film, the 
phase 6{u.) is constant throughout the surface and the 
corresponding energy vanishes. The free energy J^s to be 
minimized is a scalar generalization of Eq. |^ : 

^,, = ^Jd^u^g"f'd^e{u)dpe{u) . (3) 

Here the set of coordinates u = {ui,U2) label points on 
the surface while ^^g is the determinant of the metric 
tensor 5^/3 • On the other hand, a constant bond angle 
9{u) is not the ground state of the liquid crystal because 
it is measured with respect to an arbitrary basis vector 
Eq,(u) with a — 1,2. Indeed, it is not possible to make 
the directions of the molecules parallel everywhere on 
a curved space; the lowest energy state is attained by 
optimally distributing the unavoidable bend and splay 
of the vectors over the whole surface. The free energy 
functional JF^, to be minimized is a vector generalization 



ofEq.d) 0: 

T, ^ ^ J d^u^g''^{d^9{u)-n^{u)){di,9{u)-ilp{u)) , 

(4) 
where r2c((u), the connection, compensates for the rota- 
tion of the 2D basis vectors Eq(u) in the direction of 
Ua- Since the curl of fla{u) is equal to G'(u) 0, the 
integrand in Eq.Q) cannot be made to vanish on a sur- 
face with non-zero Gaussian curvature (riQ(u) is a non- 
conservative field and hence cannot be equal to daO ev- 
erywhere). As the substrate becomes more curved, the 
energy cost of this geometric frustration can be lowered 
by generating defects in the ground state even in the ab- 
sence of topological constraints [13, uM ■ 

In this letter, we introduce a novel coupling between a 
defect and the varying curvature of the substrate which 
originates in a conformal anomaly of the free energies of 
Equations 10) and Q. This anomaly arises, even at zero 
temperature, from imposing a constant cutoff, a, local- 
ized at the core of each defect fl9j. Bv contrast, finite 
temperature conformal anomalies |l2| are generated by 
the presence of a short wavelength cutoff for the fluctu- 
ations in 9(u) at every point on the surface. A physical 
consequence of the anomalous coupling discussed in this 
paper is that topological defects in superfluids and su- 
perconductors interact with the curvature in a radically 
different way from the case of liquid crystal order |2j . 

For thin layers of superfluids and superconductors, we 
prove that the geometric interaction -E"'(uj) is given by: 



E-'{n.)^-^qfV{u,), 



(5) 



where u^ and qi are respectively the position and topolog- 
ical charge of the defect. The geometric potential ^(u) 
satisfies a covariant version of Poisson's equation where 
the negative of the Gaussian curvature G{u) plays the 
role of the charge density: 



Dc,D"V{u) = G(u) 



(6) 



For an azimuthally symmetric surface such as the bump 
represented in Fig. ^ we can explicitly obtain V{u), as a 
function of the radial distance from the top, by employ- 
ing a two dimensional analogue of Gauss' law [111. The 
resulting potential well V{r) vanishes at infinity and its 
width and depth are given respectively by the linear size 
of the bump and its aspect ratio squared. Eq.|(Bl) has an 
elegant geometrical interpretation if a set of coordinates 
is chosen so that the metric tensor is cast in the form 
gap = Sa/3 exp {—2uj{u)) S| . The conformal factor w (u) 
is controlled by the overall shape of the surface and it 
satisfies the same Poisson Eq.® as the geometric poten- 
tial Q. We therefore proceed with the identification of 
V{u) with w(u) [23. This observation will be the basis of 
our proof of Eq.JSJ which results in the novel prediction 
that a vortex in a superfluid or superconducting film is 
repelled (attracted) by positive (negative) Gaussian cur- 
vature irrespective of its charge and sign. 




FIG. 1: A corrugated substrate and its downward projec- 
tion on a flat plane. The shaded strip surrounding P is more 
stretched than the one surrounding Q despite their projec- 
tions onto the plane having the same area. The energy stored 
in the field will be lower if the core of the defect is located at 
Q rather than P. 



For liquid crystals, the geometric interaction E'"{ui) 
contains an additional term discussed in previous inves- 
tigations of hexatic membranes 'l3|, which arises from 
the geometric frustration of the vector field. This term, 
linear in q, happens to contain the same function ^(u) 
as the conformal anomaly of Eq.©. When both con- 
tributions are included, E^{ui) acquires an unexpected 
dependence on the charge of the defect: 



E^{u,)^Kq, 



-S"'- 



(7) 



The interpretation of V^(u) as a geometric potential and 
the linear dependence on q in the first term of Eq. iQ are 
consistent with the general belief that a defect interacts 
logarithmically with the Gaussian curvature, as an elec- 
trostatic particle would with a background charge distri- 
bution. However, E^{ui) does not grow linearly with the 
charge of the defect, as expected from the electrostatic 
analogy. Instead, the geometric interaction peaks for a 
defect of charge 2tt and eventually becomes repulsive for 
q greater than the critical charge qc = 47r. Furthermore, 
E'"{ui) is symmetric around q = 2n and does not simply 
flip sign under charge inversion |22| . 

The quadratic coupling has an intuitive explanation in 
the case of azimuthally symmetric surfaces. Consider a 
very thin superfluid film deposited on the surface illus- 
trated in Fig. n with a vortex of charge q placed on top 
of the bump. In order to calculate the energy stored in 
the field, we only need to know that the superfluid phase 
0{u) changes uniformly by q along a circumference of 
length 27rr centered on the defect. Inspection of Eq.© 
reveals that the energy density of the field in the shaded 
strip at distance r is proportional to (^) , where r is 
the distance to the singularity measured in the plane of 
projection (see Fig. ^|. By vertically stretching the sur- 
face, the amount of area in the shaded strip is increased 
with respect to its projection on the plane, while the en- 




FIG. 2: Disclinations of charge 2ti located on top of bumps 
with different aspect ratios. The amount of splay in the liquid 
crystal director on the taller bump is reduced and hence the 
energy density is lower. 




FIG. 3: (a) An isolated disclination on a deformed plane feels 
a force that depends only on the enclosed Gaussian curva- 
ture, (b) The deformed plane is compactified to the sphere 
by placing a shadow defect at the south pole. 



ergy density is unchanged. As a result, the total energy 
stored in the field is greater when a vortex sits on top of 
a bumpy surface than when the same vortex is located 
at the center of a flat disk of the same area. Hence, it 
is energetically favorable for the vortex to migrate to the 
flat portions of the surface. In this case, the vortex is far 
away from the bump so that the total energy stored in 
the field does not differ much from the flat plane result 
[lj|. For less symmetric surfaces, the resulting geometric 
interaction will depend on the shape of the entire surface 
as embedded in the metric tensor. 

The physical origin of the linear coupling between de- 
fects and curvature in Eq.© is illustrated in Fig. [21 for 
a disclination of charge 2-k centered on a bump. As the 
curvature of the bump is increased, the bend or splay of 
the director of the liquid crystal decreases and hence the 
energy stored in the vector fleld is reduced. As a result, 
this linear coupling causes positively (negatively) charged 
defects to be attracted by positive (negative) Gaussian 
curvature USl. However, this mechanism competes with 
the repulsive geometric interaction illustrated in Fig. ^ 
that is at work also in the case of liquid crystal order. 
We note that the linear coupling is absent for superfluids 
because, in Eq.Q, daO{u.) is not coupled to a curvature 
dependent connection, rjQ,(u), as it is in Eq.Q. 

The critical value qc = 47r, where the single defect po- 
tential E^{\ii) changes sign, can be determined from sim- 
ple geometrical arguments. Consider an isolated discli- 
nation of charge g on a hemispherical cup placed on a 
flat plane. On account of azimuthal symmetry, the force 
acting on the defect depends only on the net Gaussian 
curvature enclosed by the circle on which it is placed, see 
Fig. |3Ji [i3 ■ This interaction is unchanged if we deform 
the outer region of the plane and eventually compactify 
it to form a sphere as illustrated in Fig. ^p. In or- 
der to satisfy topological constraints |2^ , we still need to 
place a shadow defect of charge (47r — q) at the south pole 
(the only position outside the circle that does not destroy 
the azimuthal symmetry of the initial problem). The 
curvature- defect interaction on the hemisphere is thus re- 
duced to the well known defect-defect interaction on the 
sphere [Ia|- The latter is proportional to q{^'K — q) and so 
is the curvature-defect interaction on the deformed plane 



of Fig. |3Ji, in agreement with Eq.(|7|). This provides ev- 
idence that a disclination of charge greater than An will 
be repelled from regions of positive curvature. 

We now present a derivation of the coupling between 
curvature and defects in helium and superconducting 
fllms that employs the method of conformal mapping, 
often adopted in electromagnetism and fluid mechanics 
to simplify the boundary of complicated planar regions. 
In this context, we use conformal mappings to relate the 
complex task of flnding the fleld energy on an arbitrarily 
deformed target surface to an equivalent problem on a 
homogeneous reference surface (see Fig. 0J. A conformal 
mapping has two equivalent deflning properties: angles 
map to equal angles, and very small flgures map to flgures 
of nearly the same shape. One can always find a confor- 
mal mapping from the target to the reference surface |8j 
such that qt — e^^"^"\gfl,, where gr and ga are the met- 
ric tensors on the target and reference surfaces respec- 
tively. The scaling factor e'^^"-' varies with the position u 
on the target surface, so that larger figures are inhomoge- 
neously distorted when they are mapped from the target 
to the reference surface. We choose the reference sur- 
faces to be undeformed and of the same topology as the 
target spaces (eg., g^ — Saf3 for a corrugated plane). As 
a result of the conformal transformation, defects on the 
target surface are mapped onto a set of "image defects" 
on the reference surface. The crucial property of the 
scalar free energy J^^ is its invariance under the rescaling 
of the metric by the conformal factor J2J| . However, the 
conformal symmetry of Jv, is broken upon introducing a 
short distance cutoff a that is necessary to prevent the 
energy from diverging in the core of the defect. Because 
of the varying scaling factor, the constant physical core 
radius a is stretched or contracted when projected on the 
reference space by an amount dependent on the position 
of the defect (see Fig. 0J. The radius of the image of 
the i* core is a^ = e"'-"'^a. It is this conformal anomaly 
that is responsible for generating the geometric interac- 
tion in Eq.©. In fact, the energy of the defects in the 
target space Et is equal to the energy of a configuration 
of defects (on the reference surface) whose core radii are 
position-dependent. This problem can be further trans- 
formed into the simpler task of flnding the energy En 




FIG. 4: Conformal mapping of the target surface T onto the 
reference space R. The continuous disks on both surfaces rep- 
resent the "physical" cores of constant radius a. The dashed 
Unes represent the position dependent images on R of the de- 
fect cores on T with variable radii a^. Note that the energy 
stored in the annuli comprised by the dashed and continuous 
lines in R must be added or subtracted to Er to obtain Et- 



for a set of interacting defects with constant core radius 
a plus an effective geometric potential that accounts for 
the variation of the core size with position. This geomet- 
ric potential can be derived with the aid of Fig. 01 If a^ is 
smaller (larger) than a, the energies stored in the annular 
regions indicated in Fig. 0]need to be added (subtracted) 
from Eji to obtain Et- To calculate this extra energy, 
we introduce a set of polar coordinates (r,0) centered on 
the i^^ defect. Near the defect of charge qi, the phase is 



given by « ^(j) ^iid the energy density is 



Kqf 



■. Upon 



integrating it over the annulus comprised between a and 
fli = e'^'-"'^a (see Fig. |21l, we obtain 



Et — Er 



Nd 2 

1=1 



(8) 



where Nd is the number of defects. The energy Er ac- 
counts for defect-defect interactions since any potential 
felt by a single defect would have to be constant because 
all points are equivalent on the reference surface (unde- 
formed sphere or plane). Recalling that uj{u.) = V{u), we 



recover the result of Eq.© with no dependence on the 
microscopic physics because the core size a drops out in 
Eq.ljHJ. In the case of liquid crystal order, the contribu- 
tion of the anomaly is simply added to the term linear in 
q as indicated in Eq. Q •ILj . Additional couplings result 
from abandoning the phase-only approximation to the 
superfluid energy or from higher order terms in a soft- 
spin Landau expansion. Generally, these energy contri- 
butions vary inversely as the area of the surface and can 
be made arbitrarily small by rescaling the surface, unlike 
the geometric potential in Eq.© 16] . 

Experiments that test our predictions can be realized 
by coating a bump with a thin layer of superfluid he- 
lium and rotating it around its axis of symmetry so that 
a single vortex forms |17J . The competition between 
the (repulsive) geometric interaction and the confining 
parabolic potential (generated by the rotation) would 
cause the equilibrium position of the vortex to shift from 
the center of the bump if its height exceeds a critical 
value. The vortex line could be detected by trapping of 
electrons on its core lls|. Other experiments may de- 
tect an inhomogeneous distribution of thermally induced 
defects resulting from the combined effect of the anoma- 
lous coupling and the dependence of their Coulomb-like 
interaction on the varying curvature. 

We have demonstrated that the interaction between 
defects and curvature in 2D XY-like models depends cru- 
cially on the nature of the underlying order parameter 
and we have shown how to explicitly derive the resulting 
geometric force from the shape of the surface. 
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